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Abstract 

The study of singularities and zeros of the generating functions of mul- 
tiplicity distributions is advocated. Some hints from well known proba- 
bility distributions and experimental data are given. The statistical me- 
chanics analogies provoke to look for a signature of phase transitions. The 
program of further experimental studies of the singularities is formulated. 

Multiplicity distributions in high energy collisions of various projectiles and 
targets possess qualitatively similar (but quantitatively different!) behaviour. 
That is why many fits by some well known probability distributions have been 
tried. The ever more sensitive characteristics such as the ratio of cumulant to 
factorial moments have been proposed |Q and have revealed new features of 
experimental data [0. Their understanding asks for further experimental and 
theoretical studies. It is proposed here to pay more attention to the structure 
of singularities and zeros of generating functions of multiplicity distributions. It 
is especially appealing in view of possible statistical analogies ^] . 

Let us define the generating function G(z) of the probability distribution P n 
by the relation 

oo 

G(z) = £(l + *) n P„. (1) 

n=0 

In what follows, we often use also the function 

$>(z)=lnG(z). (2) 

The (normalized) factorial (F q ) and cumulant (K q ) moments of the distribution 
P n are related to them by the formulae 

oo q 

G(z)=J2^)"F q , (F Q = F 1 = 1), (3) 



v 7 q 

$(z)^-(n)%, (#1 = 1), (4) 

9=1 

where (n) is the average multiplicity. 

First we consider some distributions which provide analytical examples for 
the nature of the singularities. We start with the fixed multiplicity (FM) dis- 
tribution when the sample of events of the same multiplicity (no) is chosen, 
then proceed to Poisson distribution (P) as a reference to independent emission 
processes and, finally, treat the gamma- (r), negative binomial (NB) and log- 
normal (L) distributions widely used to fit experimental data at high energies. 
The corresponding functions $(z) look like 

$ FM (z) = n ln(l + z), (5) 

$ P (z) = z(n), (6) 

$ r (z) = - M ln(l-Ml n (i + z )) j ( 7 ) 
M 

^ B (z) = -fcln(l-^M), (8) 

where \i and k are the adjustable parameters. The lognormal distribution is 
here the only one which is not determined by its moments. From the integral 
representation of its generating function 

<P L (z) -> - In J™ exp[- (ln * ~ ^- + x ln(l + z)]d(ln x) (9) 

it is easily seen that its convergence radius is given by the inequality 

|* + 1U<1, (10) 

i.e. the singularities come close to the point z = but they are "soft" in the 
sense that the normalization condition G(0) = 1 persists. For other distributions 
the non-trivial (essential for our purposes) singularities are situated at 

znb = k/(n), (11) 

z T = exp(/i/(n» - 1, (12) 

z P = oo, (13) 

zfm = -1- (14) 

Let us note that NB and T-singularities are close to z = if the parameters 
k and /z are much less than (n). It is especially interesting because factorial 
and cumulant moments are calculated as g-th derivatives of G(z) and &(z) at 
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that point and the nearby singularity influences their behaviour substantially. 
In particular, it is important for the ratio of the moments 

H q = K q /F q , (15) 

which is identically equal to zero for Poisson distribution, alternates sign at 
each rank in case of fixed multiplicity, and is always positive for T and NB 
tending at asymptotically large ranks to zero as q~ k ||. The different type of 
behaviour is predicted in QCD with strong decrease at low ranks followed by 
(quasi)oscillations at larger ranks jjl M, 0. It would be interesting to guess what 
singularity governs such a shape. There is no solution of the problem yet. 

Let us get some guides from experimental data. In experiments with differ- 
ent projectiles and targets the adjustable parameters are different and energy 
dependent. Nevertheless, one can get qualitative estimate of the approximate 
locations of the singular points. In e + e~ -collisions, the NB-estimates give rise 
to k/(n) ~ 1 (see, e.g., Q) and, therefore, the singularity is situated at z ee ~ 1, 
i.e. rather far from z = 0. It is much closer to the origin in /i/i-collisions where 
(see, e.g., ||) k/{n) ~ 10 _1 . The ^1^4-data is not so definite ||] (even though 
the lower statistics is slightly compensated by larger multiplicity ) and give rise 
to k/(n) < 1CP 1 and, thus, to ever closer (to the origin) singularity. The sin- 
gularities move to the origin with energy increase. Probably, these qualitative 
tendencies are related to somewhat similar regularities in the behaviour of the 
depth of the minimum of H q found for various reactions (see [|[ ^]) and to oscil- 
lations of H q at large q (see below). Moreover, the oscillations of experimental 
distributions about the smooth NB-fit (see, e.g., (t|) could be connected with 
those oscillations. Their physical meaning could correspond to various number 
of subjets (ladders etc.) contributing at different multiplicities and should be 
checked in Monte-Carlo models. Another possible source of oscillations due to 
the cut-off of the multiplicity tail by conservation laws should die out asymp- 
totically @. 

However, this cut-off plays an important role when one tries to restore the 
generating function directly from experimental data. Actually, the series (1) is 
replaced now by the polynomial in z 

N 

G N (z) = ^2(l + z) n P n (16) 

with N equal to the highest observed multiplicity. Therefore Gjv(z) has N 
complex conjugate zeros 

N 

G N (z) = l[(l-^). (17) 

j=l Z 3 

It was shown by DeWolf |ll| that the zeros cover a circle in the complex z- 
plane for ee-events generated by JETSET Monte-Carlo program at 1000 GeV. 
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It reminds of Lee- Yang zeros p2| in statistical mechanics. They do not seem 
always to close in onto the singularity of G(z) on the real axis when N tends 
to infinity. It would be interesting to check the interrelation of zeros of Gn(z) 
and singularities of G(z). 

The cumulants are determined by the moments of zeros locations 

_ (g-l)! y- (g-l)l v^ coBgflj 

\ / 3 -=i \ / 3 -=i j 

where we denote Zj — rj exp(i6j). Thus, the oscillations mentioned above are 
related to the phases of zeros. 

The study of singularities of the generating function becomes more fruitful 
if one uses statistical mechanics analogies and recalls the Feynman fluid model 
[ p"3[ H [| . The generating function is analogous to the partition function of the 
grand canonical ensemble and $(z) to free energy. The total rapidity range 
plays a role of the volume and the variable 1 + z is just the fugacity. One can 
define the "pressure" p(z) and the mean number of particles at given fugacity 
(n(z)} (proportional to the usual pressure and density) by the formulae 

p(z)=lim^, (19) 

W Z ))„ = (l + #, (20) 
az 

where §n( z ) = hiGjv(z) an d ( n (0))A r = ( n )- Let us note that the behaviour 
of (n(z))]y in the complex z-plane determined from experimental data should 
easily reveal zeros Zj of the function Gn |§ since it has poles exactly at the 
same loci Zj 

N 1 +7 

W z ))^ = E^- ( 21 ) 

3=1 

The plots olp(z) for experimental data about ee and hh- reactions extrapolated 
to Y — > oo have been shown in We have checked that the latest LEP 
data (e.g., M) well coincide with extrapolation used in ji) before the LEP data 
became available. The authors of Q claim that there is no phase transition in 
ee-collisions. The qualitative conclusion from Figs. 3a and 3b of [Q is that p(z) 
increases at z > much faster in non-diffractive ftft-collisions as compared to 
ee-collisions. It demonstrates that the /i/i-singularity is closer to the origin that 
corresponds to above conclusions. The increase would be even more drastic in 
case of AA-collisions (the data of EMU01 M were used for estimates) but it 
is strongly influenced by single events with very high multiplicity. Thus AA- 
analysis is hard to extend to large z. Probably, it has a physical origin since 
AA-collisions are the most suspected ones for phase transitions. Somewhat 
suspicious looks the constancy of p(z) at z < for /i/i-collisions in Fig. 3b of H. 
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In statistical mechanics it would be a signature for phase transition. If supported 
by further studies, it would provide hints for theoretical speculations. Really, the 
problem of phase transition in systems with relatively small number of particles 
should be treated carefully. In particular, it depends on the steepness of increase 
of p(z) with z. Some criteria of it are awaited for. However, the similarities may 
well happen to be mainly of formal nature and just the methods of analysis are 
comparable. Nevertheless, some physical models based on the analogy have 
been attempted JlJ, [if, ||, 0. 

Our preliminary qualitative results allow to formulate the further program 
of analysis of experimental data which consists of determining 

1. the radius of convergence of Gjv (1) according to Cauchi (P„ n ) and 
D'Alambert {P n /P n -i) criteria, 

2. the approach to the Carleman condition Xi^Li ^2n = oo at high ener- 
gies (N — > oo), 



3. zeros of Gn(z) (using the formulae ( |17[ ) or (|2l|)), 

4. the behaviour of the "pressure" p(z) (19), 

5. the behaviour of the "multiplicity" {n(z)) (p0[), 

6. the higher derivatives of $jv (the fractional derivatives can be used also 



18 , especially, in connection with the classification of the phase transi- 



tions of non- integer order proposed recently ||19[ ) . 

The extrapolations to Y — » oo should be attempted. It is quite probable that 
zeros locations will differ for different classes of processes (diffractive and non- 
diffractive; two- and three-jets etc). The drastic change in the behaviour of <&at 
or its derivatives must be carefully analysed to look for a possible signature of the 
phase transition. In parallel, the theoretical criteria of it in finite systems should 
be developed. We hope that the first stage of the program formulated above can 
provide some new insights into the physics of multiparticle production. More 
detailed results of it will be published elsewhere. 
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